General mechanism for amplitude death in coupled systems 
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We introduce a general mechanism for amplitude death in coupled synchronizable dynamical 
systems. It is known that when two systems are coupled directly, they can synchronize under 
suitable conditions. When an indirect feedback coupling through an environment or an external 
system is introduced in them, it is found to induce a tendency for anti-synchronization. We show 
that, for sufficient strengths, these two competing effects can lead to amplitude death. We provide 
a general stability analysis that gives the threshold values for onset of amplitude death. We study 
in detail the nature of the transition to death in several specific cases and find that the transitions 
can be of two types - continuous and discontinuous. By choosing a variety of dynamics for example, 
periodic, chaotic, hyper chaotic, and time-delay systems, we illustrate that this mechanism is quite 
general and works for different types of direct coupling, such as diffusive, replacement, and synaptic 
couplings and for different damped dynamics of the environment. 

PACS numbers: 05.45.Gg; 05,45.Pq; 05.45.Xt 



1. INTRODUCTION 

The dynamics of many complex systems can be under- 
stood as the collective behavior of a large number of dy- 
namical units coupled via their mutual interactions. The 
dynamics of such connected systems has been an interest- 
ing topic of study especially due to its relevance in under- 
standing a large variety of natural systems. Based on the 
the nature of interactions among the coupled units, they 
can exhibit many emergent phenomena such as synchro- 
nization, hysteresis-phase locking, amplitude death, and 
oscillator death [ij-yj. Among these, the phenomenon of 
synchronization is the most widely studied and has rele- 
vance in many contexts, such as neuronal networks, com- 
munication, laser systems, etc. So also, the quench- 
ing or suppression of dynamics called amplitude death is 
another emergent phenomenon of equal relevance in such 
systems. This can lead to interesting self-adjustable con- 
trol mechanisms and plays a prominent role as an efficient 
regulator of the dynamics. The occurrence of amplitude 
death has been reported in many cases, such as chemical 
reactions GruL biological oscillators |8l-ll0{. coupled laser 
systems [IJ, [l2| and relativistic magnetrons [lj|]- We 
would like to project the importance of the phenomenon 
of amplitude death in coupled systems in two contexts: 
one, as a desirable control mechanism in cases such as 
coupled lasers where it leads to stabilization [l4|, HH 
and two, as a pathological case of oscillation suppres- 
sion or disruption in cases like neuronal disorders such as 
Alzheimer's disease, Parkinson's disease, etc. [l6l - [l8j . 

The mechanisms so far reported to induce amplitude 
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death in coupled systems are de-tuning of oscillators un- 
der stron g co upling [l9l - |2lj |. coupling through conjugate 
variables |22l . |23| , dynamic coupling [24[ , and delay in 
coupling due to finite propagation or information pro- 
cessing speeds (25l-l28|. Distributed delays rather than 
discrete or constant delays have been proposed as more 
realistic models in ecology and neurobiolo gy, where the 
variance of the delay plays a relevant role [23]. So also, 
amplitude death has been studied in the context of at- 
tractive and repulsive couplings in two chaotic Lorenz 
systems [3(|- In all these mechanisms, death occurs dy- 
namically due to the targeting of the units to one or more 
of the equilibrium states or due to the stabilization of one 
of these states. The equilibrium states or fixed points can 
be either that of the uncoupled system or those evolved 
by coupling. While these mechanisms can model the am- 
plitude death observed in coupled systems of oscillators, 
we find that all these methods are system specific and 
may not work in a general case. In the case of death 
by delay coupling, the limitations of the method have 
been reported in several cases such as periodic [24| and 
chaotic systems [3ll.l32l]. Moreover, there are many cases 
such as neuronal disorders where depression of activity or 
death is due to the presence of another agency or medium 
[HI, [l7|. For such cases, the mechanism of death is still 
not fully understood and none of the above mentioned 
mechanisms so far reported is applicable. 

In this work, we introduce a mechanism for amplitude 
death caused by an indirect feedback coupling through a 
dynamic environment, in addition to direct coupling. Wc 
essentially project the role of the environment in control- 
ling the dynamics of connected systems. We find that, 
while it essentially explains quenching of activity or sup- 
pression induced by an external medium or agent, this 
method can also serve as a general mechanism to induce 
death in coupled synchronizable systems. Its generality 
lies in the fact that it seems to work in any coupled sys- 
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tern that can synchronize. It is effective in quenching 
dynamics in a variety of systems such as periodic oscil- 
lators, chaotic systems, hyperchaotic systems, and delay 
systems. We show that this method also induces am- 
plitude death in systems with different forms of direct 
coupling interactions, like diffusive, replacement, synap- 
tic coupling, etc. Specifically, we demonstrate that the 
present method, with a varied model for the environment, 
works in the case of hyperchaotic systems for which de- 
lay coupling is not effective to induce amplitude death. 
As such, it is an important step in methodology toward 
achieving controls or stabilization to desirable perfor- 
mance in many practical cases. The relevance of this 
method lies in the fact that death can be engineered and 
can be easily implemented in any system with coupled 
synchronizablc units. 

In the present work we use the indirect feedback cou- 
pling through the environment of our earlier work where 
we showed that such a coupling can induce anti-phase- 
(or anti) synchronization in two systems which are not 
directly connected [33| . Consider two systems coupled di- 
rectly such that with adequate strength of coupling they 
can exhibit synchronous behavior. Then if we introduce 
an additional indirect feedback coupling through the en- 
vironment or another external system such that it in- 
duces a tendency for anti-synchronization, then for suffi- 
cient strengths, these two competing tendencies can lead 
to amplitude death. We find that, in the state of am- 
plitude death, the subsystems stabilize to a fixed point 
of the coupled system. We also show that the method 
introduced here can induce amplitude death in coupled 
systems with different types of dynamics for the environ- 
ment and for different types of direct coupling. 

We develop an approximate stability analysis which 
provides the threshold or critical values of the coupling 
strength for amplitude death in the general context. Di- 
rect numerical simulation giving the regions of amplitude 
death in the space of coupling strengths agrees well with 
the transition curves obtained from the stability analysis. 

We also analyze in detail the nature of the transition 
to the amplitude death state. We find that all the spe- 
cific cases studied exhibit either continuous or discon- 
tinuous transitions to death. In the continuous case, as 
illustrated by two coupled Rossler systems, during the 
transition the full reverse period-doubling scenario is ob- 
served, and the system reaches a one-cycle state before 
amplitude death occurs. The transition to death then 
occurs due to a super critical Hopf bifurcation. In the 
discontinuous case, the transition is sudden due to the 
disappearance of a distant attractor and stabilization of 
a fixed point. For two coupled Lorenz systems, we find 
that the transition to death is probably via a sub-critical 
Hopf bifurcation with long transients, and prior to this, 
the systems go through a state of frustration between 
synchronized and anti-synchronized behavior. 



2. AMPLITUDE DEATH VIA DIRECT AND 
INDIRECT COUPLING 

We start with two systems coupled mutually with two 
types of coupling, namely a direct diffusive coupling and 
an indirect coupling through an environment. The dy- 
namics can be written as 

±i = f{xx) +e d /3(x2 -xi) +£elV, (la) 
X2 = f(x 2 ) + £dP{x\ - x 2 ) + e e jy, ( lb ) 
ij = -ny - y7 T (n + x 2 ). (lc) 

Here, x\ and x 2 represent two m-dimcnsional oscillators 
whose intrinsic dynamics is given by f(xi) and f{x 2 ) re- 
spectively. The systems are mutually coupled using dif- 
fusive coupling [ the second term in Eqs. ([Taj and (|lbp). 
The environment is modeled by a one-dimensional over- 
damped oscillator y with a damping parameter k. The 
environment is kept active by feedback from both systems 
as given by the last term in Eq. (fTc)) . Both systems also 
get feedback from y [ the last term in Eqs. (fTaland lTb)) ]. 
/? is a matrix (m x m) with elements and 1 and defines 
the components of x\ and x 2 that take part in the diffu- 
sive coupling. For simplicity, we take j3 to be diagonal, 
(3 = diag(P\, j3 2 , f3 m ), and in numerical simulations 
only one component f3\ is assumed to be nonzero. 7 is a 
column matrix (m X 1), with elements zero or I, and it 
decides the components of x\ and x 2 that gets feedback 
from the environment. 7 T is the transpose of 7 and it 
decides the components of xi and x 2 that give feedback 
to the environment. We take Ed to be the strength of 
direct diffusive coupling between the systems, and e e the 
strength of feedback coupling between the systems and 
the environment. 

The direct coupling Ed , gives a synchronizing tendency 
between the two systems while the coupling through the 
environment e e , gives an anti-synchronizing tendency. 
Thus, when both the couplings are above their critical 
values, there is a competition between the two tenden- 
cies, and the net result is the amplitude death. In the 
amplitude death state, the coupled systems are driven to 
a fixed point. 

We illustrate the above scheme for two coupled chaotic 
Rossler systems represented by the following equations 
(i, 3 = 1,2, i^j): 

xn = -x i2 - x i3 + £d(xji - xu) + e e y, 

Xi2 = xn + ax i2 , 

x l3 = b + x i3 (xu - c), 

y = -ny- y ^2x tl . (2) 

i 

The resulting time series for a synchronized state with 
only direct coupling, an anti-phase-synchronized state 
with only indirect coupling, and the amplitude death 
state with both direct and indirect couplings are shown 
in Fig. [TJ When e e = 0, and Ed is sufficiently large, we 



3 



X 


25 


A f\ 

WW 


WW 


x~ 


-25 


v v \J 


v/ v v v 


c\j 
X 


25 






x~ 


-25 






c\j 
X 


25 


\ w — 


(c) 1 - 


x~ 


-25 


1 1 



25 
Time 



FIG. 1: Time series of the first variables xn,i = 1,2 
of two coupled Rossler systems. (a) Synchronization for 
(£d, £e) = (0.2, 0.0). (b) Anti-phase synchronization (Ed,£e) — 
(0.0, 1.0). (c) Amplitude death for (e d ,s e ) = (0.2,1.0). Here, 
the Rossler parameters are a = b — 0.1, c = 18. The damping 
parameter of the environment is taken to be n — 1. 
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FIG. 2: Time series of the first variables xn of the coupled 
Lorenz systems for the parameters (a = 10, r = 28, b = 8/3) 
and k = 1. (a) Synchronization for (Ed,£e) = (5,0). (b) Anti- 
synchronization for (ed,£e) = (0,12). (c) Amplitude death 
for (e d ,e e ) = (5,12). 



observe synchronization [Fig.[TJa)]. When s e is increased 
for Ed = 0, the systems are in an anti-phase synchronized 
state [Fig. H^b)]. When both e e and £d are sufficiently 
large, the systems stabilize to a state of amplitude death 
[Fig. He)]. 

We apply the same scheme to two coupled chaotic 
Lorenz systems as given by the following equations: 

±n = cr{x a -Xix) +s d (xj 1 -xn) +s e y, 
x%2 = (r - x i3 )xa - x i2 , 
ii3 = xnXi2 - bx i3 , 

y = -m- Y^2 Xil - ( 3 ) 

i 

We find that amplitude death occurs in this case also. 
This is illustrated in Fig. [5J where time series for a 
synchronized state [Fig. [Ha)], anti-synchronized state 
[Fig. |2Jb)] and an amplitude death state [Fig. [He)] are 
shown. 

So far we have presented the method for identical sys- 
tems. However, the method also works for nonidentical 
systems. In general, for nonidentical systems the direct 
coupling will give a generalized synchronization between 
the coupled systems. Similarly, the anti-synchronization 
due to the indirect coupling will also become of a gen- 
eralized type. The combination of direct and indirect 
coupling still leads to amplitude death. As an example 
consider two coupled Rossler systems. In Eq. ([2]), we 
keep the parameters of one system fixed and vary the 
parameter c of the other system. We find that for suffi- 
cient strength of coupling the systems go to the ampli- 
tude death state even for large deviations in c. Also, the 
amplitude death state occurs when the individual non- 
interacting systems arc in different dynamical regimes. 
This is shown in Fig. O where time series for a generalized 



iA/WWWWWWW 


m 


"A''A'"'A''iWl' 


ft !\ i b \ i 

" /I/! L'V 

If if V i V U vi 
* 'l r*'i A A r 1 r 


i 


(c) '- 


1 1- 



60 120 

Time 



FIG. 3: Time series of the first variables xa,i — 1, 2 of two 
coupled non-identical Rossler systems. Here, the parameters 
a and b of the two Rossler systems are chosen to be same 
(a = b = 0.1), while the parameters c of the systems are 
chosen such that one of the systems is in a chaotic state (c = 
18) and the other in a periodic state (c = 4). The damping 
parameter of the environment is chosen to be k = 1, (a) 
Generalized synchronization for (ed, s e ) = (1.0, 0.0). (b) Anti- 
phase-synchronization for (ed,£e) — (0.0, 1.0). (c) Amplitude 
death for (e d ,s e ) = (1.0, 1.0). 

synchronized state [Fig. GDJa)], anti-phase synchronized 
state [Fig. [3{b)] and an amplitude death state [Fig. EHc)] 
are shown for two non-identical Rossler systems. 



3. LINEAR STABILITY ANALYSIS 

We present an analysis of the stability of the steady 
state of two systems coupled via the scheme given in 
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Eq. ([T]). For this, we write the variational equations 
formed by linearizing Eq. ([1} as 



6 = f'(x 2 )&+e d m-Z2)+e el z, 

'-6), 



~2 7 & 



(4) 



where £i, £2 and z are small deviations from the re- 
spective values. We denote synchronizing and anti- 
synchronizing tendencies through the variables £ s and £ a 
respectively as given by 



£.s = 6. - £2, 

= a + 6. 

Then Eq. (|4| can be written as 



fWlMf , f(xi)-f(*2) . 



(6) 



For stability, all the Lyapunov exponents obtained from 
Eq. ([5]) should be negative. 

In general, it is not easy to analyze the stability of the 
synchronized state from Eq. (|6]). However, considerable 
progress can be made if we assume that the time average 
values of f'(xi) and /'(X2) are approximately the same 
and can be replaced by an effective constant value fi. In 
this approximation we treat £1 and £2 as scalars. This 
approximation simplifies the problem such that only the 
relevant features remain and is expected to give features 
near the transition. This type of approximation was used 
in Refs. [H, [HJ and it was noted that it describes the 
overall features of the phase diagram reasonably well. 
Thus, Eq. (g]) becomes 
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fl^s - 2ed5,s, 

fi£, a + 2e e z, 

-K3 - y?a- 



(7a) 
(7b) 
(7c) 



We note that Eqs. (|7b|) and ([7c]) are coupled while 
Eq. (|7a[) is independent of the other two. The synchroniz- 
ing tendency is given by Eq. (|7ap and the corresponding 
Lyapunov exponent is 



(8) 



Ai = fi - 2e a 



The anti-synchronizing tendency is given by Eqs. (|7b 
and (|7c|) . The corresponding Jacobian is 



J 



and the eigenvalues are 



-e e /2 -k 



^2,3 



(M - k) ± VG" - n) 2 - 4(e 2 e - flit) 



(9) 



As noted in the preceding section, amplitude death is ob- 
tained when both synchronizing and anti-synchronizing 
tendencies are present and the corresponding coupling 
constants are greater than the critical values required 
for the respective phenomena. The synchronizing and 
anti-synchronizing tendencies become effective when the 
corresponding Lyapunov exponents, i.e., the real parts of 
the eigenvalues, are negative. From Eq. (|5|) we obtain 
the condition 

e d > fi/2, (10) 

while from Eq. (|9]) we get the following conditions: 
(1) If {fi — k) 2 < A(e 2 — fin), A2.3 are complex and the 



(5) condition of stability is 



k > fl, 



(11) 



(2) If (fi — K) 2 > 4(e^—fiK), A2,3 are real and the stability 
condition becomes 



k > fi and e e > fin. 



(12) 



If Eqs. (|T0|) and (|TT]) or (JT2|) are simultaneously satisfied, 
the oscillations can not occur and the systems stabilize 
to a steady state of amplitude death. For a given k and 
fi, the transition to amplitude death occurs at critical 
coupling strengths edc and e ec that are independent of 
each other. That is, 



and 



e dc = const, 



const. 



(13) 



(14) 



These general stability criteria are numerically verified 
for different systems in the following section. 

We can also analyze the stability of amplitude death by 
noting that the amplitude death corresponds to a fixed 
point of the coupled system. Thus, the condition for the 
stability of amplitude death is that all Lyapunov expo- 
nents of the fixed point are negative. This can be done 
for different systems numerically and is discussed in the 
next section. 



4. NUMERICAL ANALYSIS 

We apply our scheme to two chaotic systems, Rossler 
and Lorcnz. 



4.1. Coupled Rossler systems 

Now, we apply the scheme of coupling introduced in 
Eq. |T]) to the case of two chaotic Rossler systems. The 
occurrence of amplitude death in this case is illustrated 
in Fig. [TJc). This is further confirmed by calculating the 
Lyapunov exponents [35| also. When the systems are in 
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FIG. 4: Transition from region of oscillations (I) to region of 
amplitude death (II) in the parameter plane e e -Ed for coupled 
Rossler systems. Numerical simulations are done with « = 1. 
The points mark the parameter values (e ec ,£dc) at which the 
transition to amplitude death occurs. Solid triangles show the 
transition to amplitude death as Ed is increased for a constant 
e e - The horizontal line formed by these triangles confirms the 
stability condition Eq. (|13[) . Similarly, the circles correspond 
to transition to the amplitude death state as e e is increased for 
a constant Ed and confirm the stability condition of Eq. (|14fl . 



the amplitude death state, all the Lyapunov exponents of 
the coupled system are found to be negative. Figure IH^b) 
shows the largest Lyapunov exponent of the coupled sys- 
tem as a function of coupling strength e e . 

We study the transition to death by identifying regions 
of amplitude death in the parameter plane of coupling 
strengths e e — £d for a chosen value of n. To characterize 
the state of amplitude death, we use an index A, defined 
as the difference between the global maximum and global 
minimum values of the time series of the system over a 
sufficiently long interval. The case where A = repre- 
sents the state of amplitude death, while A ^ indicates 
oscillatory dynamics. The parameter value at which A 
becomes ~ is thus identified as the threshold for on- 
set of stability of amplitude death states. Using this in- 
dex, the transition curves in the parameter plane e e — Ed 
are plotted in Fig. H) We note that the points obtained 
from numerical simulations agree with the stability crite- 
ria Eqs. (|13p and (fT4|) obtained in the preceding section. 

We also verify numerically the criteria for transition 
to amplitude death given in Eq. (fT^j) . For this, the nu- 
merically obtained values of e\ c are plotted against k in 
Fig. [5] The line corresponds to the stability condition 
Eq. (|12[) and the points are obtained from numerical sim- 
ulations. It is seen that the agreement is good for larger 
values of n. However, for small values of k, the points 



FIG. 5: Transition from region of oscillations (I) to region 
of amplitude death (II) in the parameter plane k — e\ for 
the coupled Rossler system. The points of amplitude death 
are obtained numerically when the index A becomes zero. 
The solid curve is a linear fit corresponding to the stability 
condition Eq. (|12|l . with the effective fi = 0.1. The deviation 
from straight line behavior for small values of k is discussed 
in the text. 

deviate from straight line behavior. The reason can be 
seen from Eq. (fTT| which gives the lower limit on k. As 
k decreases, the damping of the environment variable y 
is reduced. However, this damping is essential for the 
anti-synchronizing tendency arising from the coupling to 
the environment. This leads to the deviations for small 
values of n. 

For two coupled Rossler systems as given in Eq. $2$, we 
study the complete phase diagram in the parameter plane 
of coupling strengths, identifying the regions of different 
dynamic states such as amplitude death, complete syn- 
chronization and anti-synchronization. Amplitude death 
states are identified using the index A as mentioned 
above. To identify synchronized or anti-synchronized 
states we use the asymptotic correlation values as the 
index, calculated using the equation 

c = < (sn(f)- < xu{t) >){x 2 i{t)- < x 21 {t) >) > 

V< (xnW 1 < xrJJ) >) 2 >< (z2i(*)- < *2i{t) >) 2 >' 

(15) 

The phase diagram thus obtained for the coupled Rossler 
system is shown in Fig. [6l When the coupling strengths 
£d and e e are small, the systems are not synchro- 
nized (white region). For small values of e e , when 
Ed is increased, the systems synchronize (light-gray re- 
gion). When e e is increased, the systems become anti- 
synchronized (dark-gray region). When both the cou- 
pling strengths are above a certain threshold as given by 
the stability conditions Eqs. (fT5|) and (HU, the systems 
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FIG. 6: Regions of different dynamical states in the parame- 
ter plane of coupling strengths E e -£d in the case of two coupled 
Rossler systems. The indices used to identify the different re- 
gions are the average correlation C [ Eq. (|15l) ] and the index 
A. The black region corresponds to the state of amplitude 
death (A ~ 0), the light-gray region to the synchronized state 
(C ~ 1), the dark-gray region to the anti-synchronized state 
(C ~ — f) and the white region correspond to the state where 
|C| < 1. Here, the parameters are the same as in Fig. [1] 



stabilize to the state of amplitude death (black region). 
We also note that the transition from complete synchro- 
nization to anti-synchronization corresponds to a phase 
transition where the average phase difference between the 
oscillators changes from to nearly w. This is similar to 
the phase- flip bifurcation reported in the context of time- 
delay coupled systems [H, H?} . 

The nature of the transitions to the state of amplitude 
death is further characterized by fixing one of the param- 
eters e e or Ed and increasing the other. This is shown in 
Fig. where the index A is plotted for increasing e e for 
a chosen value of Ed- Here, the transition from the os- 
cillatory to the amplitude death state is continuous such 
that, as the coupling strength is increased, the amplitude 
of the oscillations gradually decreases to zero. A similar 
transition is observed for the case where e e is kept fixed 
and Ed is increased. 

We also notice from the time series and phase space 
plot that, as the coupling strength increases (e<j or e e ), 
the Rossler systems undergo the full reverse period- 
doubling sequence to the one-cycle state before going 
to the amplitude death state. Then, the transition to 
the state of amplitude death occurs via a supercritical 
Hopf bifurcation. The bifurcation diagram for this tran- 
sition is shown in Fig. [8] This is further confirmed by 
computation of the fixed points of the coupled system 
and their stability near the transition region. Numeri- 
cal simulations show that the coupled Rossler systems in 
Eq. @ stabilize to the steady state corresponding to syn- 
chronized states of the subsystems. These synchronized 
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FIG. 7: The index A as a function of e e for a fixed value of 
Ed = 0.2 for two coupled Rossler systems. As s e is increased, 
we observe a continuous transition to the state of amplitude 
death (e ec ~ 0.45). 




FIG. 8: Bifurcation diagram obtained by plotting the max- 
ima of x'ii (x') for sufficiently long period by increasing the 
coupling strength e e for Ed = 0.2 and k = 1. Here, the sys- 
tems go through a reverse period-doubling bifurcation as e e 
increases, leading to a one-cycle (at e e ~ 0.35) before ampli- 
tude death occurs (e ec ~ 0.45). 



steady states are obtained from Eq. © as 



x il 


= (ci^c 2 


~Aabn/(n~E 2 e a))/2 


* 

x i2 


= -Xii/ a, 




* 
X i3 


= -VWi- 


c), 


y* 







(16) 

Of the two fixed points, the one with the plus sign in 
the second term of the x* x equation is unstable and the 
one with the minus sign in the second term of x*^ equa- 
tion become stable in the amplitude death state. The 
nature of the transition to the stable fixed point is de- 
termined by the eigenvalues of the corresponding Jaco- 
bian, and we find that at the transition, real parts of the 
complex conjugate pairs of eigenvalues become negative 



FIG. 9: (a) Real parts of the largest eigenvalues of the Jaco- 
bian obtained from Eq. , evaluated at the fixed point given 
in Eq. (|16l) . for increasing e e and fixed values of Ed = 0.2 and 
k = 1. At the transition (e ec ~ 0.435), the real parts of one 
pair of complex conjugate eigenvalues cross zero, (b) Largest 
Lyapunov exponent of the coupled Rossler system given in 
Eq. @ for increasing e e for fixed values of Sd = 0.2 and 
k = 1. The zero crossing of the largest Lyapunov exponent 
(e e c ~ 0.435) indicates the transition to the amplitude death 
state. In both figures, zero is shown as a dotted line. 



[Fig. Ufa)] , indicating a supercritical Hopf bifurcation as 
described in Ref.[38|]. In the amplitude death region, all 
the Lyapunov exponents of the coupled system [ given in 
Eq. [5] are found to be negative. The largest Lyapunov 
exponent of the coupled system crosses zero at the tran- 
sition, and this is shown in Fig. [9jb) . The nature of the 
transition is found to be the same when e e is kept fixed 
and Ed is increased. 

The above numerical results are presented for one set 
of parameters of the Rossler system. We have varied the 
parameters and verified that the method works for other 
values of the parameters. 



4.2. Coupled Lorenz systems 

We repeat the same study in the case of two coupled 
Lorenz systems. It is interesting to note that, in this 
case, the coupled systems stabilize to a fixed point that 
corresponds to anti-synchronized states for the subsys- 
tems (xn = —X2i,x\2 = -x 2 2,xi3 = Z23) as shown 
earlier in Fig. [He). The regions of different dynamical 
states in the parameter plane of coupling strengths in 
this case are shown in Fig. [TO] When both Ed and e e are 
small, the systems are not synchronized (white region). 
For very small values of Ed and large E e , the systems are 
anti-synchronized (dark-gray), and when Ed is increased 
from this state, the systems go to the amplitude death 
state (black region). For small values of £ e and large 
Ed, the systems are synchronized (light- gray ) . As s e in- 
creases, the systems first lose synchronization, and for 
larger values of e e , they stabilize to the state of amplitude 
death (black). In the de-synchronized state before the 
amplitude death state, the attractor in the phase space 
is highly distorted and the system goes through a state 
of frustration, trying to stabilize to the anti-synchronized 
state from the synchronized state before death occurs. 
This is illustrated in Fig. [TTJ where the time series of the 
synchronization error between the two Lorenz systems is 




FIG. 10: Regions of different dynamical states in the param- 
eter plane of coupling strengths e e Sd in the case of Lorenz 
systems. The black region corresponds to the state of ampli- 
tude death (A ~ 0), the light-gray region corresponds to the 
synchronized state (C ~ 1), the dark-gray region correspond 
to the anti-synchronized state (C ~ —1) and the white region 
correspond to the state where |C| < 1. 
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FIG. 11: Synchronization error (iu — X21) as a function of 
time in the case of Lorenz systems for Ed = 5 and e e = 4.665. 



shown. Near this transition region, some initial condi- 
tions remain in a chaotic transient state for a long time 
before becoming stabilized to the fixed point. The phe- 
nomena of multi-stability and hysteresis are also observed 
in this region. The nature of the transition to amplitude 
death in this case is shown in Fig. [12J Unlike the case of 
coupled Rossler systems, here we see that the amplitude 
of oscillations drops suddenly at a critical strength of cou- 
pling. Thus, the transition is directly from the chaotic to 
the amplitude death state. A similar type of transition 
from the chaotic to the amplitude death state in the case 
of time-delay coupled Lorenz systems has been reported 
in Ref. [26|. We further characterize this transition by 
computing the fixed points of the coupled system given 
in Eq. ([3]) and evaluating their stability near the transi- 
tion region. 

Numerical simulations show that the coupled Lorenz 
systems in Eq. ((3J) stabilize to the steady state corre- 
sponding to the anti-synchronized state of the subsys- 
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FIG. 12: The index A as a function of e e for a fixed value of 
Ed = 5.0 for two coupled Lorenz systems. Note that the tran- 
sition to amplitude death is sudden (e ec ~ 7.88) as opposed 
the case of coupled Rossler systems where the transition is 
continuous. 



tems. These steady states are obtained from Eq. Q as 



' ii 



' 12 



= ±1 



((r - l)a - 2s d )b 



(T + 2e d 



(a + 2e d ) 



Xl3 


ab 


* 


* 


X 21 






* 


X 22 


— ~~ x Y2i 






X 23 




* 

V 


= 0. 



(17) 

For both the solutions, we find that at the transition, 
the real parts of the complex conjugate pairs of eigen- 
values of the corresponding Jacobian become negative 
[Fig. 113( a)]. For e d > £ dc , numerically an unstable limit 
cycle is found to coexist with the stable state of ampli- 
tude death for certain initial values. As there is no stable 
limit cycle before amplitude death, and an unstable fixed 
point becomes stable, it seems that this is a sub-critical 
Hopf bifurcation. All the Lyapunov exponents of the 
system [Eq. [3] are found to be negative at the amplitude 
death state. At the transition, the largest Lyapunov ex- 
ponent of the coupled system becomes negative as shown 
in Fig. [□]» . 





FIG. 13: (a) Real parts of the largest eigenvalues of the 
Jacobian obtained from Eq. Q, evaluated at the fixed point 
given in Eq. (|17|) . for increasing e e and fixed values of Ed — 5 
and k = 1. At the transition (e ec ~ 7.31), the real parts of one 
pair of complex conjugate eigenvalues cross zero, (b) Largest 
Lyapunov exponent of the coupled Lorenz system given in 
Eq. ([3])for increasing E e for fixed values of Ed — 5 and n — 1. 
The zero crossing of the largest Lyapunov exponent (e ec ~ 
7.43) indicates transition to the amplitude death state. In 
both figures, zero is shown as a dotted line. 



a variety of systems and find that, for all cases which are 
synchronizable, the method works. We understand that 
extensive numerical simulations may not analytically es- 
tablish the generality of the mechanism. However, the 
results of our numerical simulations appears to indicate 
that this method is quite general. 

We also give an intuitive physical argument to support 
our claim of generality of the method. Our mechanism 
consists of having two types of coupling. The first is di- 
rect coupling, which leads to synchronization. We know 
that, if the coupling constant is sufficiently large, the syn- 
chronization condition ensures that the largest Lyapunov 
exponent transverse to the synchronization manifold is 
negative. Considering the space of coupled oscillators 
as a product of the individual systems and a network 
of two nodes, the synchronization manifold corresponds 
to the direction e s = (1, l) T in the network coordinates. 
Similarly, the coupling through the environment which 
ensures anti-synchronization leads to the condition that 
the largest Lyapunov exponent transverse to the direc- 
tion e a = (1,— 1) T is negative. Since we have coupled 
only two systems, ensuring that the largest exponents 
transverse to both synchronizing and anti-synchronizing 
directions, i.e. e s and e a , are negative, implies that all 
the Lyapunov exponents are negative. Thus the system 
must converge to a fixed point. 

In this section, we present the results of applying this 
method to periodic, time-delay, hyperchaotic, and driven 
systems and different schemes of direct coupling. 



5. AMPLITUDE DEATH IN OTHER CASES 

We have presented a mechanism for inducing ampli- 
tude death in coupled systems due to the competing ef- 
fects of synchronizing and anti-synchronizing tendencies. 
We have shown this in the context of two chaotic sys- 
tems, namely, coupled Rossler and coupled Lorenz sys- 
tems. To test the generality of the method, we apply it to 



5.1. Amplitude death in periodic systems 

We study two standard limit cycle oscillators, namely 
Landau-Stuart and van der Pol oscillators, coupled using 
the scheme given in Eq. ([I]). The Landau-Stuart sys- 
tem is a nonlinear limit cycle oscillator, which has been 
previously used as a model system for studying the phe- 
nomenon of amplitude death [22|, [H|. In our case, the 
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dynamics of two coupled Landau-Stuart systems is given 
by the following set of equations: 

± a = [x\ x + x%)x a - ujx i2 + e d (xj! - x a ) + e e y, 

Xi2 = {xfi +X- 2 )Xi2 +CJXU, 

y = -Ky- y^2 x ^- ( 18 ) 

i 

From numerical analysis of the above equations with 
uj = 2, we see that, for small values of e e and Ed, the sys- 
tems are synchronized. For small values of Ed and large 
values of s e , the systems are in the anti-synchronized 
state. When the strengths of both Ed and e e are suffi- 
ciently large, we observe amplitude death. We note that 
here the transition to the state of amplitude death is sud- 
den, similar to the case of Lorenz systems. However, in 
this case, for a given strength of coupling, the stability of 
amplitude death state depends on the initial conditions 
indicating multi-stability, that is, some initial conditions 
go to the amplitude death state, while some other initial 
conditions remain in the oscillatory state. Such a multi- 
stability has also been reported for amplitude death phe- 
nomena in the case of Landau-Stuart oscillators using 
conjugate coupling [22 |. For some initial values, ampli- 
tude death occurs even in the absence of direct coupling 
( Ed = 0) ■ A possible explanation is that the fi of the 
individual system is negative or zero such that, the stabil- 
ity condition k > [i given in Eq. pip is always satisfied. 
It is seen that the area of the basin of amplitude death 
increases as Ed is increased. 

We repeat the same study for the case of two coupled 
periodic van der Pol systems given by the following equa- 
tions: 

in = x i2 + £d{x 3 i - xa) + e e y, 
xi2 = a(l - x 2 ll )x i 2 - Xn, 

y = -Ky (19) 



5.2. Amplitude death in time delay systems 

The Mackcy- Glass time-delay system is well studied 
as a model exhibiting hyperchaos. Stabilization to fixed 
point, or amplitude death in such systems has been re- 
ported by use of stabilization methods such as conven- 
tional feedback, tracking filters and delayed feedback [3!| . 
Here we consider two Mackey- Glass systems coupled via 
both direct and indirect couplings as given by 

Xi = -axi + ^ Tl m + £d(x 3 - Xi) + £ e y, 

- 1 - ' X ri 

y = -*y- y x j> ( 20 ) 

.7=1,2 

where Xi represents the Mackey-Glass system [4£j and x T i 
is the value of the variable Xi at a delayed time Xi (t — r) . 
The parameters of the Mackey-Glass systems are chosen 
such that, the individual systems are in the hyperchaotic 
regime ( a = 1, /3 = 2, r = 2.5, m = 10). For very small 
values of Ed and e e , the two systems are not synchronized. 
The systems synchronize as Ed is increased. For large 
values of e d , if E e is increased, the systems go to a state 
of amplitude death. When Ed is small and e e is large, 
the systems are in an anti-synchronized state. As we in- 
crease Ed, the systems go to the amplitude death state. 
However, for large values of Ed and s e , the Mackey-Glass 
systems appears to show a different behavior from the 
Rossler or Lorenz systems. We observe a re-entrant be- 
havior to rhythmogenesis, both as e e increases and as Ed 
increases. This transition also satisfies our stability con- 
ditions Eqs. (TT5)) and (|bfj) . We find that the transition to 
the state of amplitude death in the case of two coupled 
Mackey-Glass systems is continuous and that the systems 
go through a reverse period-doubling sequence reaching 
a limit cycle before the amplitude death occurs. This is 
similar to the case of Rossler systems discussed earlier. 



We choose the parameter a = 1 such that the sys- 
tem has a stable limit cycle when both the couplings are 
absent (i.e., e e = Ed = 0). For a suitable strength of 
the direct coupling Ed, the systems synchronize and am- 
plitude death is observed when both direct and indirect 
couplings are above their respective thresholds. As far 
as the nature of the transition to the state of amplitude 
death is concerned, the van der Pol oscillator shows a 
different behavior from that of other systems. We find 
that the nature of the transition to amplitude death de- 
pends on the type of coupling parameter. Fixing Ed and 
increasing E e , we see a smooth transition similar to that 
in the case of Rossler systems and by fixing e e and in- 
creasing Ed, we get a sudden transition, as in the case of 
Lorenz systems. 



5.3. Amplitude death with alternate schemes for 
direct coupling 

So far, we have studied direct coupling of the diffusive 
type. Synchronization is also possible with direct cou- 
pling of different types. Here, we now study two such 
types of coupling. 



5.3.1. Lorenz systems with replacement coupling 

Here we consider Lorenz systems coupled using a dif- 
ferent scheme of coupling, namely replacement coupling, 
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as given by the following equations 



V 



a(x j2 - xn) +e e y, 
(r - xa)xii - Xi2, 

X t lX t 2 — bXi3, 



-Ky 



2 



E 



xn. 



(21) 



Here, the direct coupling is of the replacement type, 
such that the X2 variable in the first function of the first 
system is that of the second system, and vice versa. This 
type of coupling leads to synchronization as reported in 
Ref. [!l||. We introduce indirect coupling through the 
variable y. We find that, for suitable values of the cou- 
pling strength e e , the systems stabilize to a state of am- 
plitude death. 



5.3.2. Synaptically coupled Hindmarsh-Rose model of 
neurons 

The Hindmarsh-Rose system is a model of neurons 
which shows the spiking and burstin g be havior of the 
membrane potential of a single neuron [42j . We take two 
neurons with excitatory synaptic coupling [43[ between 
them and introduce an indirect coupling as given by the 
following equations: 



Xi2 
Xi3 

y 



%i2 - x iX + ax a - x l3 + I + 



V r - xn 



1 + exp(-A(ar 3 -i - 9)) ' 



1 - bx a - x i2 , 
p(s{xii +x) -x i3 ), 



-Ky 



tE 



xn. 



(22) 



0.5 




FIG. 14: Region of amplitude death (black) in the parameter 
plane of coupling strengths e e -£d for two coupled Hindmarsh- 
Rose neurons. We choose the parameters of the system to be 
a = 3,6 = 5, p = 0.006, s = 4, x = 1-6, and I = 3.2. The 
parameters in the coupling terms are chosen to be V r = 3, 
6 = —0.25, A = 10, and k = 1. For the neuronal systems, the 
criterion used in the text for identifying the amplitude death 
state, i.e., the index A ~ 0, needs to be changed. Here, we use 
the maximum voltage from both neurons (x m ) as the index 
to identify the death state. The criteria for synchronized and 
anti-synchronized states also needs modification. Hence, the 
other regions are not shown explicitly in the figure. 



is possible with other intrinsic dynamics for the environ- 
ment also. 

For this, we consider the case of two Rossler systems 
coupled with a dynamic environment, where the intrinsic 
dynamics of the environment is that of a damped har- 
monic oscillator. The dynamics of the coupled system in 
this case is given by 



Here, the variable xn represents the membrane poten- 
tial of a neuron and the variables x\2 and xa are related 
to ion currents across the membrane. We choose the pa- 
rameters of the system such that the individual neurons 
are in the chaotic bursting state. When e e = 0, and 
is sufficiently large, the bursts of both neurons become 
synchronized. For larger values of Ed and e e , we observe 
patches of amplitude death (shown in Fig. H4)) . Thus as 
we keep one of the coupling parameters fixed (e e or Ed), 
and increase the other, we observe a transition to the 
amplitude death state and again re-entrant behavior to 
spikes. 



5.4. Amplitude death with alternate dynamics for 
the environment 

In the previous sections, we have taken the intrinsic dy- 
namics of the environment to be that of an overdamped 
harmonic oscillator. Here, we show that amplitude death 



xn = —Xi2 - x i3 + £d(xji - xn) + e e yi, 

Xi 2 = xn + axi2, 

x i3 = b + x l3 (xn -c), 

Vl = V2 - -^^xn, 

i 

m = -yi - Ky x . (23) 

Here, the variables y\ and y2 represent a two-dimensional 
environment together forming an underdamped harmonic 
oscillator. For very weak coupling (Ed ~ 0,e e ~ 0), 
the two Rossler systems are not synchronized. When 
the coupling strength Ed is increased while e e is kept 
fixed at zero, the systems become synchronized. On the 
other hand, when e e is increased while Ed is kept fixed 
at zero, the systems become anti-phase-synchronized. 
When both £ e and Ed are above their respective thresh- 
olds, amplitude death is observed. 

We repeat the same study by taking the intrinsic dy- 
namics of the environment as that of an overdamped 



11 



Duffing oscillator. The equations in this case are 

111 = -Xi2-x i3 + £d(xj 1 -x i i)+e e y, 

112 = xn + ax l2 , 

x l3 = b + x l3 {x l i -c), 

i 

For small values of e e and Ed, we see that the systems 
are not synchronized. They become synchronized as 
Ed is increased from this state. For small values of Ed 
and large values of e e , the systems are in an anti-phase- 
synchronized state. When the strengths of both Ed and E e 
are sufficiently large, we observe amplitude death. The 
phase diagram in this case is qualitatively similar to that 
given in Fig. [S] 

5.5. Amplitude death in hyperchaotic Rossler 
systems 

We also consider the case of two hyperchaotic Rossler 
systems as given by the following equations 

Xil = -Xi2 - x i3 + £ e yi, 

— Ed cos 0(sin 0(xii — xji) + cos9(xi 3 — 2^3)), 
X12 = xn+ ax l2 + xa + £ e y 2 , 

X i3 = b + X i3 Xii 

—£dsin6(sin6(xii — Xji) + cos9(xi 3 — Xj 3 )), 
Xi4 = -cx i3 + <rx i4 + E e y 3 , 

Vi = -nyi - y^an-i, 

i 

yi = -Ky 2 - y ^2,x l2 , 

i 

m = -Ky 3 - y ^x l4 , (25) 

i 

where i,j = 1,2 and j 7^ i. 

We choose the parameters of the system such that, the 
intrinsic dynamics of the systems is hyperchaotic. For 
this system, the method of time-delay coupling is found 
ineffective for producing amplitude dcath[32[. In the ab- 
sence of coupling via the environment, the direct coupling 
via a scalar signal results in synchronization of the two 
hyperchaotic systems for suitable values of parameters Ed 
and 9, as reported in Ref. [44|. We take the environment 
to be three-dimensional in this case. We find that am- 
plitude death occurs in the coupled system for suitable 
values of the coupling strengths £ e and n. This is shown 
in Fig. [HI 

5.6. Small oscillations in driven systems 

We apply the scheme described in this paper to driven 
systems such as driven van der Pol and Duffing systems. 
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FIG. 15: Time series of the first variables xn,i = 1, 2 of the 
two coupled hyperchaotic Rossler systems given in Eq. (|25p 
in the amplitude death state. Parameters of the systems are 
chosen to be a — 0.25, 6 = 3, c = 0.5, and a = 0.05. The 
parameters used in the direct coupling term are Ed = 2.5, and 
9 — 7r/3 and in the indirect coupling terms are n = 10 and 

Ee =4. 

In such driven systems, the fixed point is not a solution 
for the individual or coupled systems. So the interpreta- 
tion of amplitude death as in other systems needs to be 
changed. Here, we interpret the amplitude death state 
as the state of very small amplitude oscillations. 

Driven van der Pol systems with direct diffusive cou- 
pling and indirect coupling through the environment can 
be written as 

Xil = Xi2 + £d(xjl - Xn) + £ e y, 

Xi2 = a(l - x a )xi2 —xn+P cos(wt), 

y = -ny- y^^i- (26) 
3 

We find that, when £ e = and Ed increases, the two 
systems become synchronized. In this state, if we start 
increasing £ e , we get a state of small oscillations or am- 
plitude death. Figure ITBT a) shows the time series for 
such a state. 

Figure [TTT a) plots the index A as a function of e e . We 
first see a transition from a limit cycle to two different 
limit cycles for the two systems. This state subsequently 
goes to the amplitude death state continuously as e e in- 
creases further. On the other hand, if we keep e e fixed 
and increase Ed, we find a sudden transition to the am- 
plitude death state. As Ed increases further, we find a 
continuous transition to a state of increasing amplitude 
oscillations. This is shown in Fig. [T7Tb). 

The same study is repeated for the case of two coupled 
Duffing systems given by the following equations: 

in = x i2 + £d(xji - xn) + e e y, 

Xi2 = -ax i2 + xn — ar?! + /3 cos(wt), 

y = -Ky- y J2 X ^- ( 2? ) 
3 

We find that a regime of small oscillations is possible in 
this case also. This is shown in Fig. [ToTR - ). 

Thus, wc have illustrated that the method for inducing 
amplitude death in coupled systems introduced in this 
paper works for periodic, time-delay, hyperchaotic, and 
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FIG. 16: Time series of the first variable xn of the coupled 
driven systems showing small oscillations. Here, we interpret 
the small oscillations as the state of amplitude death (see 
text), (a) Driven van der Pol system for (ed,£e) = (1.0,3.5). 
The parameters of the individual systems are taken to be 
a = 8.53, cj = 0.63, and j3 — 1.2 (b) Driven Duffing system 
for (sd,£ e ) = (1.0,4.0). The parameters of the individual 
systems are taken to be a = 0.25, u = 1, and /? = 0.3 



driven systems. It is effective in quenching the dynam- 
ics even with different forms of direct coupling, such as 
replacement coupling and synaptic coupling of neurons. 

The context of coupled neurons presents an important 
case of amplitude death which could explain the mech- 
anism of disruption or suppression of synaptic signals 
in the case of neuronal disorders, such as Alzheimer's 
disease, as being due to induced activity and feedback 
through a protein called amyloid beta (A/3). Using nu- 
merical studies on a few neuronal models, we have shown 
that the competing effects of synaptic activity and the 
indirect interaction mediated by the protein A/3 lead to 
sub-threshold activity and synaptic silencing |45| . 



6. DISCUSSION 

In this paper, we show that indirect coupling through 
a dynamic environment in addition to direct coupling 
can lead to amplitude death in chaotic systems such as 
Rossler and Lorenz systems. The approximate stability 
analysis developed for general cases gives the transition 
region in parameter space which is further supported by 




FIG. 17: Transition to amplitude death state in two coupled 
driven van der Pol oscillators. Here the death state corre- 
sponds to a state of small oscillations since the fixed point is 
not a stable state. The index A remains finite although very 
small in the amplitude death state, (a) The index A as a 
function of e e for fixed Ed = 1.0. The transition is continuous 
(Eec ~ 3.08). There is a sudden transition due to a pitchfork 
bifurcation of the limit cycle observed at e e ~ 2.55. (b) The 
index A as a function of Ed for fixed e e = 3.5. The transition 
to the amplitude death state is sudden (sdc ~ 0.72). We also 
observe a re-entrant continuous transition to periodic oscilla- 
tions with increasing amplitude at Ed ~ 1.14. 



direct numerical simulations. The nature of the tran- 
sition to amplitude death is found to be of two typical 
types, one continuous and the other discontinuous. 

In conclusion, the method for determining amplitude 
death introduced in this paper is quite general and 
works for different types of systems such as periodic, 
chaotic, hyperchaotic and time-delay systems and also 
with different types of direct coupling, for example, dif- 
fusive, replacement coupling, coupling via scalar signals 
and synaptic coupling. We have demonstrated that our 
method works in hyperchaotic Rossler systems where 
time-delay coupling is ineffective in inducing amplitude 
death. We have also presented a physical argument for 
the generality of the method. Although the theory is 
developed for the case of identical systems, we have veri- 
fied that this method of inducing amplitude death works 
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in the case of non- identical systems as well. In fact, we 
did not find any exception to our scheme provided the 
coupled systems are synchronizable. 

We also note that the method introduced in this paper 
can be implemented in practical cases. What is needed is 
the design of a suitable environment which can introduce 
the appropriate indirect coupling between the systems. 
Moreover, in many natural systems, the environment or 
external medium exists and can be instrumental in caus- 
ing suppression of the dynamics. The mechanism pre- 



sented here provides an explanation for this phenomenon 
in such cases. 
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